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ABSTRACT 



We compute the temperature and IR signal of particles of radius a and 
albedo a at heliocentric distance R, taking into account the emissivity effect, 
and give an interpolating formula for the result. We compare with analyses 
of COBE DIRBE data by others (including recent detection of the cosmic IR 
background) for various values of heliocentric distance, R, particle radius, a, 
and particle albedo, a. We then apply these results to a recently-developed 
picture of the Kuiper belt as a two-sector disk with a nearby, low-density sector 
(40 < R < 50 — 90 AU) and a more distant sector with a higher density. We 
consider the case in which passage through a molecular cloud essentially cleans 
the Solar System of dust. We apply a simple model of dust production by comet 
collisions and removal by the Poynting-Robertson effect to find limits on total 
and dust masses in the near and far sectors as a function of time since such 
a passage. Finally we compare Kuiper belt IR spectra for various parameter 
values. Results of this work include: (1) numerical limits on Kuiper belt dust as 
a function of (R, a, a) on the basis of 4 alternative sets of constraints including 



those following from recent discovery of the cosmic IR background by [Hausei 
et al. (1998)| ; (2) application to the two-sector Kuiper belt model finding mass 
limits and spectrum shape for different values of relevant parameters including 
dependence on time elapsed since last passage through a molecular cloud cleared 
the outer Solar System of dust; and (3) potential use of spectral information to 
determine time since last passage of the Sun through a giant molecular cloud. 
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1. Introduction 



Two papers examining IR flux upper limits from COBE DIRBE and IRAS ( |Backma~ 



et al. 1995| ; |Stern 1996a| ) have put quite stringent limits on the total amount of mass in 



the form of dust (particles less than about one centimeter) in the Kuiper belt (KB) region 
inside ~ 50 AU. Roughly, they find that the dust mass, M D , is bounded by M D < 10" 5 M e . 
In both cases, however, the results are based on specific, albeit quite reasonable, models 
for the distribution of dust as a function of particle radius, a, and heliocentric distance, R, 
and an assumption as to albedo, a. The objectives of this paper are: (i) to address KB IR 
radiation from an arbitrary distribution; (ii) to obtain COBE DIRBE limits on dust mass, 
Mo, as a function of a, R, a; (iii) to apply these results to a case of particular interest for 
a variety of parameter values; and (iv) to take into account the recent analysis of COBE 
DIRBE data by [Hauser et al. (1998)| in which bounds at 60 and 100 microns, and values at 
140 and 240 microns, for the cosmic infrared background have been found (see Table 1). 
We note that recent reviews of the Kuiper belt include Stern (1996c) and Weissman and 
Levison (1997). 

The case of interest explored here is that in which the near portion of the Kuiper belt 
(40 < R < [50 — 90] AU) is significantly depopulated. Such a structure has been shown 
( [Stern 1996b| ) to make it possible to understand formation of the large bodies observed 
( Pewit t fc Luu 19"§5| ; |Jewitt, Luu fc Chen 1996] ) in the face of the large eccentricities 



observed, which would imply erosive rather than growth-promoting collisions in the present 
day 30 — 50 AU region. The explanation ( [Stern 1996c| ) appears to be an originally more 



massive region of low eccentricities to permit growth of larger bodies in a short accretion 
time, with the region later suffering physical and dynamical erosion, and depopulation due 
to dynamical interaction with Neptune ( Levison fc Duncan 1993 ; Holman fc Wisdon 1993 ; 



Malhotra 1995; [Duncan et al. 1995| ; |Morbidelli et al. 19961 ; Levison et al. 1997). In this 
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picture, only the component of the Kuiper belt sufficiently far past these effects would 
retain the surface density distribution, R~ 2 to _R _3//2 , of the original solar nebula. The 
near belt (e.g., 30 — 70 AU) would be significantly depopulated. A key implication of this 
result would be different IR signals from each of the two sectors since, on the average, the 
dust in the far sector would be cooler, dimmer, and in orbits of much smaller eccentricity. 
Accretion in the early Kuiper belt is also studied by |Kenyon fc Luu (1998) . 

The plan of the paper is as follows. In Section 2 below, we give our results on dust mass 
limits as a function of heliocentric distance, particle radius, and particle albedo (R, a, a) 
from COBE DIRBE IR signal limits. In Section 3 we apply these results to the mass limits 
for the two sector distribution discussed above. In Section 4 we discuss the variation of 
spectra with model parameters. Section 5 contains our concluding discussion. 



2. IR Dust Signal 

The IR intensity signal in the wavelength band A to A + dX from iV particles, each of 
radius a and temperature T at heliocentric distance R, spread uniformly in an ecliptic band 
from —9 to +8 is given by 

I(N, 6, a, R, T, X)dX = (na 2 / R 2 ) [N/(4tt sin9)]e IR I BB dX (1) 

here I BB is the radiation emitted by a black body per unit area-time- wavelength-solid angle: 

I BB (X, T) = 2hc 2 /[X 5 (e hc / XkT - 1)] (2) 

the emissivity factor, em, which gives the suppression of radiation of wavelength larger than 
the particle is assumed to be: 

e IR (a,X) = [l,(a/X) n ] [X<a,X>a]. (3) 
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We investigate below the implications of choosing n to be 1 and 2. A useful discussion of 
emissivity and also absorptivity factors is given by Packman fc Paresce [1993] . Because IR 



observations integrate their view radially (i.e., view the radial distribution in projection on 
the sky), in applying results based on equation ([!]) to specific models, one must integrate 
over the model dependent distribution in radius a and heliocentric distance R. We will do 
this in Section 3 below. The temperature for each particle can be found from its values of 
a, R and a. Equating solar power absorbed with power emitted gives the standard result 

tt a 2 (l - a)L /(47r R 2 ) = 4tt a 2 / /(A, T)e IR (X)dX (4) 

Jo 

where L is the luminosity of the Sun. We introduce a variable y by 

T = hc/(kay). (5) 

Equation (f|) can then be written in the form 

z = (1 - a)L Q a A /(?,2Ti 2 hc 2 R 2 ) = F(y) (6) 

Below we invert equation (§) to find y, and hence the temperature which is needed to use 
equation (JJ) to find the IR signal from the dust as a function of z. In equation (H) we have 
(for n=l) 

F(y) = [yh(y) + h(v)]/y 5 (7) 
where the functions fk are given by 

/■oo 

fk(y)= x k - l dx/[e x -l}. (8) 
J y 

Inverting equation @ can be done analytically for z near zero and infinity, with the results 

y=(f A (0)/z) 1/4 (9) 

near infinity, and 

y = (/ 5 (oo)/z)^ (10) 
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near z = 0. In equations (||) and (10), the asymptotic forms are given in terms of the 



ordinary gamma function, T, and the Riemann zeta function, (, by 

f k = r(k)c(k). (n) 

We have evaluated the integrals in equation @ numerically and found for n=l that a simple 
quintic approximation to Iny represents y to better than one percent over the interval [] 
15 > Inz > —15. The approximation formula is 

y = e Ao z B (12) 

where 

B = J2Mlnz) n (13) 
i 

with A = 0.4624; A 1 = -0.2433; A 2 = -0.001846; A 3 = 0.0001243; A A = 2.799 x 10~ 6 ; 
and A$ = —3.066 x 10~ 7 . A similar procedure can be used for other values of n in equation 
(H). |Backman fc Paresce (1993~)| give a different expansion, but the temperature predictions 



are the same. 

With the results of equations (Q, [FT], |T^) for particle temperature, we can use equation 
(ffl) to find the total IR signal from a given mass M in the form of particles of radius a, 
density p, and albedo a located at heliocentric distance R. We insert in equation ([]]) that 

N = M/(4vrpa 3 /3). (14) 



We can then compare the results of equation ([!]) with COBE DIRBE data. In Section 
3 we will integrate over a distribution in R and a in the model considered. Here we provide 
"model-independent," unintegrated limits . We make comparisons with "data" four times: 
first we choose, as the constraint on the KB IR signal, the full COBE DIRBE signal (as 



4 For higher and lower z equations (|9] and [It]) are preferred. 
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given in Figure 2 of [Backman et al. 1995|) ; second, we choose the constraints that result from 
(mo del- dependent) subtraction of the signal from the asteroid belt as cited by Backman et 
al. (credited to, e.g., Reach 1992, 1988). Our third choice is the COBE DIRBE two sigma 
upper limits for A =60 and 100 microns and values for 140 and 240 microns found by |Hauser 



et al. ( 1998)| . Finally we choose the set of standard deviations from [Hauser et al. (1998 



Our treatment of the Hauser et al. results is motivated by the fact that their procedure for 
subtracting the Solar System contribution to the IR signal did not include explicitly fitting 
and subtracting a Kuiper belt contribution, so that the implications of their work for the 
KB are not clear. Further discussion of their treatment of Solar System source subtractions 
is given in Pwek et al. [1998] and |Kelsall et al. ( 1998 )| . The values in the four bands for 



each of the four sets of constraints are given in Table 1. 

In Figures 1 and 2 we display the results of using equation (|l|) to find the maximum 
mass M(a, R, a) (in units of M©) of material in the form of particles of radius a with 
albedo a at heliocentric distance R. We use a particle density of 1 gm cm -3 and take 
the KB density distribution as constant in density to ± 1/3 radians and zero outside this 
region. Our results can, of course, be scaled to other values for these parameters. Figure 1 
gives results for a = 0.05, which would be the case if the dust is either made mainly from 
the surface of comets or made sufficiently long ago that it has suffered the same radiation 
darkening as is believed responsible for the low albedos of comet nuclei (e.g. |Johnson 1990| ). 



For comparison we also give results in Figure 2 for a = 0.5, which might obtain if the dust 
were made of relatively pure ices. 

In Figure la we give the curves for all four constraint sets of Table 1 for a = 0.1 cm. 
These curves do not depend on the value of n in equation (3) because little radiation is 
emitted of wavelength larger than 0.1cm. In Figure lb we give the results for particle 
radius a = 0.001 for all four constraint sets and for n = 1 and n = 2. Figure 2 gives the 
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cor responding curves for a = 0.5 cm. Finally, Table 2 gives dust temperatures, for n = 1 
and n = 2 for selected values of a, R, and a. 

Note that: (1) Figures 1 and 2 give the total mass that, if located exactly at heliocentric 
distance R (and optically thin), would give an IR signal that would exhaust the constraint, 
not a distribution or a mass density; (2) the bound on the dust mass is linear in the assumed 
bound on the signal whereas the dust mass itself and the IR signal from it, in models for 
dust production such as those of Section 3 below, are proportional to the square of the total 
mass of colliding comets that gives rise to the dust (as noted in Packman et al. 1995|) ; (3) 
our model for dust is unsophisticated compared to work addressing dust properties in more 
detail such as Liou et al. (1996) , Ishimoto (1998)| , and Yamamoto fc Mukai 1998] ; (4) in the 



absence of the emissivity effect of equation (|3|), the mass limit M(a, R, a ) scales linearly 
with a (the ratio of volume to surface area). The emissivity effect can raise or lower mass 
limits according to the relative values of long and short wavelength signal limits, since its 
effect is to shift energy from long wavelength bands to short; and (5) the places in Figures 
1 and 2 at which there are discontinuities in slope are places at which there is a change in 
which wavelength band is limiting. 



Some of us, in other work flAnderson et al. 19981 ) have looked at Kuiper belt information 



that can be derived from Pioneer 10 which has been in the region of the KB for over a 
decade. A limit on the mass of the KB in the form of objects of radius somewhat under a 
centimeter of roughly a tenth of an Earth mass (averaged over the region of the KB out 
to about 65 AU) has been found based on the fact that Pioneer encountered no object 
energetic enough to penetrate its propellant tank. That result is clearly compatible with 
the results of this paper (i.e., that result is also an upper bound although one less stringent 
than the IR bound of this paper). 
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3. Two Sector Belt Model 

We now apply the results above to the Kuiper belt model outlined in the Introduction, 
in which there is an under populated near region and a denser far region. We assume 
that the boundary is somewhere between 50 and 90 AU depending on such factors as the 
influence of Neptune on orbits and the effects of velocity evolution in the disk. For this 
paper we attempt to bracket the boundary by presenting results for it being at 50 AU and 
at 90 AU; where we need a representative value, we choose 70 AU. 

We insert the results of Section 2 for the IR signal as a function of heliocentric distance 
R and particle radius a into a simple computer model for dust production. Our model for 
such production is similar to those described previously flBackman et al. 1995) ; |Stern 1996a| ), 



but extends them to take into account dependence on time since passage through a GMC 
last cleansed the Solar System of dust. We assume that the region is populated, at first, 
solely by planetesimals (comets) of 5 km radius, compute the collision rate, assume that all 
collisions result in a dust distribution with an a -3 5 dependence, (in Section 4 we consider 
a -3,0 as well) and that each dust particle is removed at its appointed Poynting-Robertson 
(P-R) time. More precisely, the model is as follows. 

The total collision (dN co i(R) / dt) rate is found from the expression Nin 2 uv for collisions 
per unit time of N\ particles with a distribution n 2 per unit volume of others, with a the 
collision cross section and v the relative velocity. For us, n 2 is equal to nc, the number 
density of comets, and N% = nc^n R 2 dR sin 9. For the collision of 5 km comets it is 

d 2 N col (R)/dRdt = n 2 c (R)7r a 2 c v c (R) (4vr R 2 sin9) (15) 

In equation (|15l), ric(R) is the number of comets per unit volume, a>c the 5 km radius, and 
vc the (average) relative velocity. We take a surface distribution, of total mass Mt between 



-10- 



Ri and R2, falling like R 2 to obtain 

n c (R) = (M T /Mc)/[ln(R 2 /R 1 )47Tsin6R 3 } (16) 

where Mc is the mass of an individual comet (about 5 x 10 17 gm for ours, assuming unit 
density) and 9 is the half-height extent of the KB in declination. We take the average 
relative velocity in equation fll5|) to be given in terms of the average eccentricity e by 



v c = V2v K {R) 



5 9 ,9 

- < e > + < 1 > 
4 



1/2 



(17) 



where vk is the Kepler velocity, the \/2 factor sums in quadrature the velocity deviations 
from circular of the two colliding objects, and the sum in square brackets is from 
Lissauer and Stewart (1993). Below, we make the standard equilibrium approximation 

< % >=< e > /2. 

In our model, the IR signal from the KB region between Ri and R2, from particles of 
radii between a\ and a 2 is given by 

I(X) = / dR da[^N M {lCi/dBM][dN D /da]T{R,a,t)h{a, R, 6, a). (18) 

J Rl J a\ 

In equation (|T8|), 1\ is the wavelength distribution of equation ([[]) for iV = 1 (i.e. for each 
particle). Here, r is the shorter of: (i) the time t since last passage through a molecular 
cloud eliminated Kuiper belt dust ( [Stern 1990|) ; and (b) the Poynting- Robertson time, tpR, 
for a particle of radius a at heliocentric distance R. (See [Burns, Lamy fc Soter 1979 | for a 



full, lucid discussion of this and related effects). The P-R time, tpR, is given by 

t PR = 47rc 2 paR 2 /{3L G (l- a)] (19) 
which may be evaluated approximately by 

t PR ^ 10 10 yrp(gmcm- 3 )a(cm)(R(AU)/A0) 2 /(l - a). (20) 
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Other effects that need, in principle, to be included are radiation pressure and erosion 
by the ISM. Radiation pressure is only important for a on the order of a few microns where 
the P-R time is already very short so little error is made in neglecting it. Stern (1990) has 
shown that most of the ISM effect comes from occasional passage of the Sun through a giant 
molecular cloud (GMC). We assume that the ISM effect can be approximated by periodic 
removal of dust from the Solar System and compute dust production as a function of time 
since the most recent such "tidying." There is an important question with regard to this 
assumption: GMC passage erodes about 10 gem -2 from comets in the Oort cloud flStcrn 



1990|) and for ones in the KB. For a KB of 0.1M® in 5 km comets this implies creating a 
dust mass of about 10 _5 M e - of the order of the maximum dust mass compatible with 
IR signal bounds (Backman et al. 1995, Stern 1996). However the high relative velocity 
between the Solar System and the GMC (about 30km/s) means that the dust generated 
by the collisions will tend to have high enough velocity (roughly 10% of the 30km/s) to 
exit the Solar System flFujiwara et al. 1989|) so that much less bound dust will be generated 
(and that which is generated will be expelled in the next encounter). 

We assume this Solar System bound dust is negligible. This assumption is conservative 
in the sense that, if some of this dust were captured into solar orbit, our upper limits (below) 
on total KB mass would be made lower (more constraining). We take the distribution 
dNjj/da to be the one that is self replicating under collisions (Dohnanyi 1969 and, more 
recently in more generality, [Tanaka et al. 1995|) ; in terms of mass it is given by m~ n / 6 , 



while in terms of radius, we have 

dN D /da = n a- 7/2 . (21) 

In Section 4, we briefly discuss the effect of varying from equation (|2T|) . We assume density 
p = lgm/cm 3 for both comets and dust and we use equation (|2lJ) up to 5 km. Then, by 
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normalizing equation fl2*T|) to the total volume of the two comets, we have (with a in cm) 

n = a 5 J 2 /2 « 2 x 10 14 . (22) 

We have coded equations (fL5|)-(|2~2|). For each set of initial parameters, we evaluated 
the IR signal in the four bands for a low value of total KB sector mass, and then increased 
the sector mass until the IR limit was exceeded in one of the four bands, 60, 100, 140, 240 
microns (for one of the four constraint sets of Table 1). The results are given in Tables 3 
and 4 and Figures 3 and 4. In the tables, the columns are: (1) "case number;" (2) signal 
constraint applied (1-4 as in Table 1); (3) albedo a; (4) n in equation (§); (5) inner radius 
of sector; (6) outer radius of sector ; (7) time since dust began accumulating; (8) maximum 
sector mass in Earth masses; (9) total dust produced (one micron to 0.5 centimeter) in units 
of 10 _5 M e ; (10) dust remaining after Poynting-Robertson effect in the same units. Table 3 
gives results for 64 sets of input parameters for the near sector of the Kuiper belt and Table 
4 for the far. In each, the first 32 cases are for n = 1 the last 32, n = 2. We have taken two 
values for the time r since the KB was last cleaned of dust: the age of the Solar System and 
200 million years. The latter value is about twice the minimum needed for comet nuclei 
radiation mantles to darken their albedos to the 0.05 average observed in the few cases for 
which there are measurements flJohnson 1990|) . As noted above, we take 50 and 90 AU as 



extreme possibilities for the boundary between the near and far sectors. As before, we take 
the KB comet number density to be constant in declination to 9 = ±1/3 radians for the 
near sector, For the far sector, which did not suffer dynamical interaction with Neptune as 
did the near, we use < e >= 9 = 0.01 radians. 

We turn first to Table 3 and the near sector. Features of the calculations and results 
that might be noted include: 

-As expected, the IR signal rises as the square of the total belt mass, Mt (or the comet 
number density), since the collision rate does. 
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-As < e > (or 6) is increased, the relative velocity goes up linearly, but the density 
goes down linearly so the total collision rate in equation fll5|) stays constant; the signal 
in equation (|18D then varies as Mf./e because of the 6 dependence of I\. Thus the table 
can be extended approximately to other eccentricity values (within our simple model) by 
M T oc (e) 1 / 2 . 

-We have good agreement with the results of Backman et al. Case 15 in Table 3 
(Mt = 1.17 for a = 0.5, t = 4.5 Gyr) corresponds to their parameter and constraint 
choices but needs to be corrected for their choice of p = 0.5, 9 = 1/6, and boundaries 
of 30-100 AU for them versus 40-90 AU for case 15. This gives a limit of about 
M e /2 [= 1. 17(1/2) (l/ v / 2)(10/9) 3 / 2 ] which is their result (Af©/3) within better than a factor 
of 2. 

-The ratio of dust remaining (after P-R loses) to total produced, M DR /M D , goes up 
with increasing R 2 because the Poynting- Robert son time goes up as R 2 . The masses Mdr 
and Md each scale as < e > since they vary as M|.. 

-There is a factor of about 30 between the largest near sector upper mass limit, 7.0 M ffi 
(case 40), and the smallest, 0.22 (case 29), but all values are above current best estimates 
of about 0.1 M e for the near sector. 

-There is a factor of about 70 between the largest calculated value of remaining dust, 
about 41 x 10" 5 M e (case 39), and the smallest, about 0.56 x 10 5 M e (case 26). 

-We have given results for a = 0.05 and a = 0.5 for comparison although it is unlikely 
that the latter can obtain in reality. 

-Note the relative insensitivity to the emissivity index, n; going from n = 1 to n = 2 
raises the temperature of small grains significantly, but it suppresses the rate of emission 
from them in the wavelength range above their size; these two effects tend to cancel. 
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-We have chosen the upper limit for particle radius at a 2 = 0.5 cm, but one can 
understand the dependence of the IR signal on a 2 simply. For n = n a _/3 , the IR signal goes 
as noa~^aHpR(a)da ~ a 2 _/3 for tpi?(a 2 ) < t the present age of the Solar System. The 
IR signal stops growing with a 2 when tpp(a 2 ) = t . From equation (0) we see that, for 
p = 1, the signal stops growing with a 2 at a 2 = 0.5 cm for a <C 1 (as noted by Backman et 
al. 1995). 



-Note that we can see, from Equations (|l5l), fll8|) , (|21|) , and (j22|) , the effect of varying 
our assumption of 5 km for the comet radius, ac- We have, by equation ( p2|) n ~ a^ 2 and 



d 2 N col /dRdt ~ M 2 ja% which gives / ~ M 2 /d% or M T 



3/2 „ „ fl 3/4 



-As emphasized by Backman (1995, 1998), there is a strong dependence of the IR flux 
on the assumption made for the inner edge of the dust distribution. In our calculations, we 
assume that it is at 40 A.U. where Duncan, Levison and Budd (1995) find the approximate 
boundary between stable and unstable orbits (Backman et al (1995) and Stern (1996) 
assume a hard inner edge at the orbit of Neptune). Were it the case that the P-R effect 
operated unimpeded, then, in the steady state, the signal from KB dust in the inner Solar 
System would swamp the signal from KB itself. Estimating the amount of dust that 
survives into the inner Solar System is difficult. Important effects include: dynamical orbit 
stability as a function of grain mass; collisions with interstellar dust (see Liou, Zook, and 
Dermott 1995); and, particularly, sublimation of volatiles even out as far as 30 A.U. Thus 
the results of this section on IR bounds on the KB mass should be interpreted as quite 
conservative. A substantial signal from KB dust in the inner Solar System would mean 
much stronger bounds on the KB mass and/or the time since last GMC encounter. 

Table 4 gives the corresponding results for the far sector of the KB. We take the outer 
boundary of the far sector to be 120 AU and, again, use 50 and 90 AU for the (inner) 
boundary with the near sector. One sees that the model limits on total sector mass are 
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similar to those on the near sector. The fact that they are not larger is a result of several 
factors including: e and 9 are assumed down by a factor of 30 (cutting My by a factor 
of v^O ); P-R dust losses fall quadratically with R; and the emissivity effect of equation 
(§) causes temperatures to fall less rapidly than R~ l l 2 for small a. It is very important 
to note that the dust production model may not apply to the far sector. This is because, 
for low values of e, one expects collisions to be adhesive rather than fragmenting in nature 
with little dust production. The definition of "low," however, depends on the unknown 
(but size- dependent strength of Kuiper materials. As noted above, the limits for < e > 
high enough to result in fragmentation can be found from the results given by using the 
< e > x l 2 scaling law. Thus, for example, if we take 100 m/s as a rough lower limit for 
strong materials, then for < e > greater than 0.033 at R = 100 AU the results of Table 4 
(and Figure 3b and 4b below) scaled by [< e > /0.01] 1 / 2 would obtain. 

We turn now to consideration of the time dependence of the IR signal. In Figure 3 
we give the upper bounds on the near and far sector masses for the four constraint sets of 
Table 1. Here we choose the sector boundary at 70 AU. Again we plot the results for the 
case of eccentricity e = 1/3 for the near sector and e = 0.01 (assuming collisions result in 
fragmentation) for the far sector. In computing the IR signal we have used, for albedo, 
a = .05 at large times and a = 0.5 at small times with the break at = 10 8 years in 
order to see maximum possible time variation. However, it is unlikely that even passage 
through a giant molecular cloud could lead to an a value as high as 0.5. In computing the 
Poynting-Robertson losses we used an interpolating function for a in equation(|19|) 

a (t) = 0.5 -0.45(1 -e" t/tD ). (23) 

We see from Figure 3a that near sector, upper mass bounds decrease by about a factor 
of 2 from short times (10 7 years) to 10 s years and about another factor of 5 to long times 
(4.5 x 10 10 years.) The bounds on the far sector, Figure 3b, fall somewhat more steeply. 
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In Figure 4 we plot total dust produced and dust remaining for the two sectors as a 
function of time since last passage through a molecular cloud for the case of constraint 
set 2. One sees from the divergence of the "total produced" and "dust remaining" curves 
that, in time, the total mass of dust removed by the P-R effect becomes considerable. In 
Section 4 we will see that P-R removal of small particles, which by the emissivity effect of 
equation (^) have higher temperatures, results in a change over time in the spectrum at 
long wavelengths. This is because the amount of high T material stays constant while the 
amount of low T material grows. 

A final effect of interest is potential destruction (Liou et al. 1996) of KB dust grains of 
radii 4.5 — 25 pm by 26km/s ISM grains in the mass region 10 -15 — 10~ 12 gm discovered by 
the space craft Ulysses ( prun et al. 1993|) . If the flux of such grains is relatively constant 



over time and position, removal of KB grains in the 4.5 — 25 pm range would decrease 
the IR signal significantly for the case t=tpR. For example, if it has been 10 8 years since 
last GMC passage, this effect could significantly decrease the signal from 80 AU or more 
(assuming that < e > is high enough to produce dust). 



4. Spectra 

In this section we explore a little, within the two sector model, the extent to which 
spectral information in a Kuiper belt signal could provide information on belt structure and 
properties or, contrariwise, on how the KB spectrum varies with belt parameters. First, 
Figure 5 shows the spectrum, of each sector, for three different values of p (0.5, 1.0, 2.0). 
The calculations are for n=l and a=0.05 and sector mass of 0.1 M®. Note that doubling 
p halves the signal for small wave lengths, but quarters it for long wavelengths. This is 
because short wavelength radiation is emitted by small particles heated by the emissivity 
effect of equation (|3]). Doubling p for fixed sector mass results in halving the number of 
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comets and hence quartering the number of collisions and the amount of dust. Since it 
also doubles the dust P-R lifetime, the signal from small particles is only halved. For large 
wavelength, the signal is quartered because it comes largely from particles with tpp already 
near, or greater than, to- 

Figures 6 (n = 1) and 7 (n = 2) give spectra for the near and far sectors, for 0.2 and 
4.5 Gyr since molecular cloud passage, and for two values of (3 in the particle distribution 
with radius (n a _/3 ). We choose, for the two values of (3, (3 = 3.5 as used in Section 3 and 
(3 = 3.0. this latter value was used by Weissman and Levison (1997) to fit observations 
and bounds for sizes over 1 km. Thus it is somewhat speculative to apply it in the sub-cm 
range. Nevertheless it provides a measure of the variation in bounds and in spectra with (3 
and it is of some special interest since it sits on the dividing line between signal domination 
by large particles and signal domination by small particles. For /3=3 the signal integral, 
J*" 2 daa'^a 2 , for a 2 > apR {tpR(app) = t ) varies logarithmically with a 2 - We use for a 2 , in 
this case, the same 0.5cm as for (3 = 3.5 for comparison purposes, rather than the 10 km 
above which the Weissman-Levison fit has (3 changed to 4.5 (the correction would be an 
overall factor of about 10). 

In Figures 6 and 7, we have chosen 70 AU for the boundary between near and far 
sectors, e = 1/3 and 0.01 respectively for the two sectors, p = 1, a mass of 0.1 M® for both 
sectors in all cases, and albedo a = 0.05. Figure 6 has emissivity index n = 1 in equation 
(0); Figure 7, n = 2. The first observation to be made with regard to the results is the 
striking variation among the spectra according to time since last passage through a GMC 
The 0.2 Gyr and 4.5 Gyr intensity curves agree for A < 20 microns, but differ by over an 
order of magnitude for A > 100 microns. Therefore the IR Kuiper belt signal spectrum 
would appear to hold significant promise for determining the time since such passage. 

We also see that the (3 = 3.0 spectra, for short wavelength, fall below the (3 = 3.5 
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spectra by about roughly the expected value (10 km/1 cm) 1 / 2 ~ 10 3 where 1 cm is the 
effective upper limit for the a integral for (3 = 1 and 10 km is the upper limit for (3 = 3. 
Note that the factor of 33 in the ratio of the values chosen for < e > for the two sectors 
roughly cancels the effect of the greater distance to the far sector. One might also recall 
from the discussion of Section 3 that the intensity curves scale with M T , < e >, and a c as 
M|/(< e > a J 2 

In Figure 8, we compare spectra for times of 10 7 , 2 x 10 8 , and 4.5 x 10 9 years for near 
and far sectors (using: (3 = 3.5, n = 1; a = 0.04; M = 0.1 M e ; and < e >= 0.333, 0.01 for 
near and far sectors respectively). Again the ratio of the values of < e > chosen makes the 
intensities for the two sectors similar to within an order of magnitude. In both cases we 
see that for wavelengths around 30 microns, the ratio of the intensities for the two extreme 
times is about a factor of 4 while that for A > 200 fim is over 100. Future measurements 
of the KB IR spectrum should provide, in addition to much compositional information, 
information on the length of time since the belt was last cleaned of dust. 

As in Section 3, we have assumed that the inner edge of the KB dust distribution is at 
40 A.U. If the P-R effect, in spite of the dust destruction effects such as sublimation cited 
there, were to result in substantial amounts of KB dust pentrating to the inner Solar System 
there would be corresponding modifications in the spectra as viewed from the Earth. 



5. Conclusion 



It is tempting, in thinking about IR constraints on the Kuiper belt, to work to the 
result ( |Backman et al. 1995| ; |Stern 1996a|) of about 10 _5 M e of dust of radii between 1/im 
and 5 mm. However the results of Sections 2 through 4 show that rather wide ranges 
of values for total Kuiper belt mass (factor of 30) and for belt dust (factor of 70) are 
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possible. Tables 3 and 4 and Figures 1-3 show that there can be major differences in this 
mass limit according to assumptions made about albedo, distribution in particle size, and 
heliocentric distance, subtraction of IR signal from non-belt sources, emissivity effect index, 
belt inclination, and time since last passage through a molecular cloud. As discussed above, 
our results are in the approximation that the effect of the ISM can be approximated (Stern 
1990) by removal of Solar System dust upon passage through a giant molecular cloud. 
Survival of dust from such an encounter and/or subsequent ISM destruction or creation 
of KB dust could change our conclusions, but we do not believe these effects will make 
major changes for the reasons given above. As noted in Section 3, our mass limits are quite 
conservative in that we do not include the IR signal from dust particles that might, under 
the P-R effect, penetrate to the inner Solar System in spite of such dust destruction effects 
as sublimation of volatiles. 

If a surface mass density of the form S(-R) = Xo-R -2 is normalized to the mass of solids 
in the outer planets between 5 and 35 AU (about 48 M e ) we have S ~ 3.9M e AU~ 2 . 
Such a distribution would give about 13M e between 40 and 70 AU and the same amount 
between 70 and 120 ( and a little over half that much between 90 and 120). One can see 
from Figure 3 that the IR limits make it unlikely that the near sector of the Kuiper belt 
can have the full Xo-R -2 surface density distribution of the inner planets. If the far sector is 
dynamically cool enough that collisions are adhesive rather than fragmenting, there is little 
dust and no real IR limit on sector mass. 

The spectra studied in Section 4 tend to indicate that future spectral information could 
help distinguish Kuiper belt parameters. In particular, one result of the Poynting-Robertson 
effect is that such spectrum characterizations as the ratio of long wavelength intensity to 
short wavelength intensity could be a good measure of time since last passage of the Sun 
through a giant molecular cloud. 
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Fig. 1. — (a) Model-independent mass limits (in units of M®) for particles of radius 
a = 0.1 cm and albedo a = 0.05 as a function of heliocentric distance R for the 4 sets of 
constraints of Table 1. The curves for the case n — 1 in equation (|3]) are essentially identical 
to those for n = 2 as explained in the text. The digit by each curve is the constraint set 
number, (b) Limits (in units of M e ) for particles of radius a = 0.001 cm and albedo a = 0.05 
as a function of heliocentric distances for the 4 sets of sector constraints. Curves with n = 1 
in equation (^) are not identical to ones with n = 2 as in Figure la because the equilibrium 
temperatures for large bodies at the heliocentric distances involved produce significant IR 
with A > a = 0.001. The pair of digits by each curve gives constraint number and n. 

Fig. 2. — (a) Limits for a = 0.1 cm with a = 0.5, as explained in the caption for Figure la. 
(b) Limits for a = 0.001 cm with a = 0.5 as explained in the caption for Figure lb. 

Fig. 3. — (a) Upper limit on near sector mass as a function of time since last passage 
through a GMC for the 4 different constraint sets of Table 1 (in units of Earth masses). 
Note < e >= 0.333. (b) Upper limit on far sector mass (as in Figure 3a). Note < e >= 0.01. 

Fig. 4. — Dust produced and dust remaining (after P-R effect) as a function of time since last 
passage through a GMC (in units of M e ) for both near and far sectors assuming saturation 
of the mass bound for constraint set 2 of Table 1. Here < e > is 0.333 for near sector and 
0.01 for far. 

Fig. 5. — Spectrum /(A) variation with density p (of both comets and dust) for near and far 
sectors. Note, as explained in text, I ~ for short A but p~ 2 for long. 

Fig. 6. — Spectrum variation with (i) sector [S=(N,F) = (near, far)], (ii) time since last 
GMC passage [t=(R,D) = (recent, distant) = (2 x 10 8 , 4.5 x 10 9 )yr], and (iii) number density 
distribution index (a _/3 ,/3 = 3.0,3.5) for albedo a = 0.05 and n = 1 in equation 3. Curves 
are labeled by (S,t,/?). For reference, in descending order at A = 0.1 cm we have: (F,D,3.5), 
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(N,D,3.5), (F,D,3.0), (F,R,3.5), (N,R,3.5), (N,D,3.0), (F,R,3.0), (N,R,3.0). 

Fig. 7. — Spectrum variation as in Figure 6 for n = 2. Note that the values at A = 0.1 cm 
are just as in Figure 6 (as might be expected since essentially all long wavelength radiation 
comes from a > A which is unaffected by the emissivity effect). 

Fig. 8. — Spectrum variation with (i) time since last GMC passage [t=10 7 , 2 x 10 8 , and 
4.5 x 10 9 yr ] and (ii) sector [S=(N,F)] for (3 = 3.5, a = 0.05, and n = 1. 
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Table 1. COBE DIRBE Upper Limits on the IR Signal From the Kuiper 

Belt in MJysr" 1 . 



A (cm) 


T a 


I b 


J c 


/ d 


0.006 


16.0 


0.3 


1.50 


0.75 


0.010 


8.0 


1.0 


1.27 


0.63 


0.014 


5.2 


2.5 


1.17 


0.32 


0.024 


3.1 


2.0 


1.12 


0.24 



''Constraints from data at —1/3 radians declination. 



b Constraints subtracting the asteroid belt contribution ( Backman ct al. 1995 ) 



Constraints from Hauser et al. (1998), CIB two -sigma limits (A=0.06 and 0.01) and 
values (A=0.014 and 0.024). 

Constraints from |Hauser et al. (1998)j CIB one sigma errors. 
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Table 2. Dust Temperature As A Function of R, a, a and n. 

R T T T T T 

(AU) a= (1 cm) (0.1cm) (0.01cm) (lO^cm) (10" 4 cni) 

n=l a=0.05 

40 43.9 43.9 45.7 67.5 106.9 

60 35.9 35.9 38.0 57.4 90.9 

80 31.1 31.1 33.4 51.1 81.0 

100 27.8 27.8 30.3 46.8 74.1 

120 25.4 25.4 28.0 43.5 68.9 

n=l a=0.5 

40 37.4 37.4 39.4 59.3 94.1 

60 30.6 30.6 32.9 50.5 80.0 

80 26.5 26.5 29.0 45.0 71.3 

100 23.7 23.7 26.4 41.1 65.2 

120 21.6 21.6 24.4 38.2 60.6 

n=2 a=0.05 

40 43.9 43.9 46.6 86.2 185.7 

60 35.9 35.9 39.1 75.3 162.2 

80 31.1 31.1 34.6 68.4 147.4 

100 27.8 27.8 31.6 63.5 136.8 

120 25.4 25.4 29.4 59.8 128.8 
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Table 2 — Continued 



R T T T T T 

(AU) o= (1 cm) (0.1cm) (0.01cm) (10- 3 cm) (lQ- 4 cm) 



n=2 a=0.5 

40 37.4 37.4 40.5 77.5 166.9 

60 30.6 30.6 34.2 67.7 145.8 

80 26.5 26.5 30.4 61.5 132.4 

100 23.7 23.7 27.8 57.1 123.0 

120 21.6 21.6 26.0 53.7 115.7 
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Table 3. Total and Dust Masses for Near Sector. 



Case 


Type a 


a 


n D 


Rl 


R'2 


t 

(Gyr) 


(Mffi) 


M D 

(io- 5 m ) 


M DR 
(10- 5 M e ) 


[ 




0.05 




40.0 


50.0 


4.5 


0.86 


55.43 


26.37 


2 




0.05 




40.0 


50.0 


0.2 


1.75 


10.18 


9.17 


Q 
O 




u.uo 






on n 
yu.u 


A £ 
^±.o 


z.oy 


fifi Qfl 
uu. yu 




4 




0.05 




40.0 


90.0 


0.2 


5.52 


13.46 


12.29 


5 




0.50 




40.0 


50.0 


4.5 


0.80 


47.28 


29.47 


g 




0.50 




40.0 


50.0 


0.2 


1.95 


12.59 


11.72 


7 




0.50 




40.0 


90.0 


4.5 


2.44 


58.94 


39.29 


8 




0.50 




40.0 


90.0 


0.2 


6.18 


16.84 


15.83 


9 


2 


0.05 




40.0 


50.0 


4.5 


0.30 


6.50 


3.09 


10 


2 


0.05 




40.0 


50.0 


0.2 


0.50 


0.81 


0.73 


11 


2 


0.05 




40.0 


90.0 


4.5 


0.95 


8.95 


4.79 


12 


2 


0.05 




40.0 


90.0 


0.2 


1.59 


1.11 


1.02 


13 


2 


0.50 




40.0 


50.0 


4.5 


0.36 


9.73 


6.07 


14 


2 


0.50 




40.0 


50.0 


0.2 


0.60 


1.18 


1.10 


15 


2 


0.50 




40.0 


90.0 


4.5 


1.17 


13.66 


9.11 


16 


2 


0.50 




40.0 


90.0 


0.2 


1.91 


1.62 


1.52 


17 


3 


0.05 




40.0 


50.0 


4.5 


0.42 


13.12 


6.24 


18 


3 


0.05 




40.0 


50.0 


0.2 


0.82 


2.21 


1.99 


19 


3 


0.05 




40.0 


90.0 


4.5 


1.28 


16.38 


8.77 


20 


3 


0.05 




40.0 


90.0 


0.2 


2.59 


2.95 


2.70 


21 


3 


0.50 




40.0 


50.0 


4.5 


0.41 


12.69 


7.91 


22 


3 


0.50 




40.0 


50.0 


0.2 


0.92 


2.83 


2.64 
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Table 3 — Continued 



Case 


Type a 


a 


n b 






t 

(Gyr) 


M be lt 

(Me) 


(io- 5 m ) 


M DR 
(10- 5 M e ) 


23 


3 


0.50 


1 


40.0 


90.0 


4.5 


1.28 


16.34 


10.89 


24 


3 


0.50 


1 


40.0 


90.0 


0.2 


2.93 


3.79 


3.56 


25 


4 


0.05 


1 


40.0 


50.0 


4.5 


0.22 


3.58 


1.70 


26 


4 


0.05 


1 


40.0 


50.0 


0.2 


0.43 


0.63 


0.56 


27 


4 


0.05 


1 


40.0 


90.0 


4.5 


0.67 


4.48 


2.40 


28 


4 


0.05 


1 


40.0 


90.0 


0.2 


1.37 


0.83 


0.76 


29 


4 


0.50 


1 


40.0 


50.0 


4.5 


0.22 


3.49 


2.17 


30 


4 


0.50 


1 


40.0 


50.0 


0.2 


0.48 


0.77 


0.72 


31 


4 


0.50 


1 


40.0 


90.0 


4.5 


0.67 


4.46 


2.98 


32 


4 


0.50 


1 


40.0 


90.0 


0.2 


1.53 


1.04 


0.97 


33 




0.05 


2 


40.0 


50.0 


4.5 


0.88 


57.71 


27.46 


34 




0.05 


2 


40.0 


50.0 


0.2 


1.88 


11.73 


10.56 


35 




0.05 


2 


40.0 


90.0 


4.5 


2.65 


69.88 


37.40 


36 




0.05 


2 


40.0 


90.0 


0.2 


5.99 


15.81 


14.43 


37 




0.50 


2 


40.0 


50.0 


4.5 


0.82 


49.75 


31.01 


38 




0.50 


2 


40.0 


50.0 


0.2 


2.19 


15.91 


14.81 


39 




0.50 


2 


40.0 


90.0 


4.5 


2.50 


62.30 


41.53 


40 




0.50 


2 


40.0 


90.0 


0.2 


7.02 


21.77 


20.47 


41 


2 


0.05 


2 


40.0 


50.0 


4.5 


0.30 


6.70 


3.19 


42 


2 


0.05 


2 


40.0 


50.0 


0.2 


0.51 


0.87 


0.79 


43 


2 


0.05 


2 


40.0 


90.0 


4.5 


0.96 


9.16 


4.90 


44 


2 


0.05 


2 


40.0 


90.0 


0.2 


1.65 


1.20 


1.09 
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Table 3 — Continued 



Case 


Type a 


a 


n b 




R 2 


t 

(Gyr) 


(M e ) 


M D 
(10- s M e ) 


(10" 5 M e ) 


45 


2 


0.50 


2 


40.0 


50.0 


4.5 


0.37 


9.93 


6.19 


46 


2 


0.50 


2 


40.0 


50.0 


0.2 


0.62 


1.27 


1.18 


47 


2 


0.50 


2 


40.0 


90.0 


4.5 


1.18 


13.87 


9.24 


48 


2 


0.50 


2 


40.0 


90.0 


0.2 


1.98 


1.73 


1.63 


49 


3 


0.05 


2 


40.0 


50.0 


4.5 


0.43 


13.60 


6.47 


50 


3 


0.05 


2 


40.0 


50.0 


0.2 


0.87 


2.48 


2.23 


51 


3 


0.05 


2 


40.0 


90.0 


4.5 


1.31 


16.98 


9.09 


52 


3 


0.05 


2 


40.0 


90.0 


0.2 


2.76 


3.37 


3.07 


53 


3 


0.50 


2 


40.0 


50.0 


4.5 


0.42 


13.23 


8.25 


54 


3 


0.50 


2 


40.0 


50.0 


0.2 


1.02 


3.48 


3.24 


55 


3 


0.50 


2 


40.0 


90.0 


4.5 


1.31 


17.12 


11.41 


56 


3 


0.50 


2 


40.0 


90.0 


0.2 


3.30 


4.79 


4.51 


57 


4 


0.05 


2 


40.0 


50.0 


4.5 


0.22 


3.73 


1.77 


58 


4 


0.05 


2 


40.0 


50.0 


0.2 


0.47 


0.72 


0.65 


59 


4 


0.05 


2 


40.0 


90.0 


4.5 


0.68 


4.64 


2.48 


60 


4 


0.05 


2 


40.0 


90.0 


0.2 


1.49 


0.97 


0.89 


61 


4 


0.50 


2 


40.0 


50.0 


4.5 


0.22 


3.63 


2.26 


62 


4 


0.50 


2 


40.0 


50.0 


0.2 


0.54 


0.98 


0.91 


63 


4 


0.50 


2 


40.0 


90.0 


4.5 


0.69 


4.68 


3.12 


64 


4 


0.50 


2 


40.0 


90.0 


0.2 


1.74 


1.34 


1.26 



a Table 1 constraint number. 
b See equation (3). 
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Table 4. Total and Dust Masses for Far Sector. 



Case Type a a n b Ri R 2 t M bdt M D M DR 













(Gyr) 


(M ffi ) 


(10- 5 M ffi ) 


(10- 5 M e ) 


1 


1 

1 


n n^ i 
U.UO J 


£n n 
L OU.U 


i on n 
IzU.U 


A K 
4.0 


n C7 
U.o / 


O OO 

z.yz 


l.oO 


z 


1 

1 


n n^ i 
U.UO J 


£n n 
L OU.U 


i on n 
IzU.U 


n o 
U.z 


o no 
z.uy 


n 7ft 
U. ( 


n 71 
U. / 1 


Q 

o 


1 

1 


n n^ i 
U.UO 


on o 
l yu.u 


i on n 
IZU.U 


A K 
4.0 


1.0 f 


7 /I 
f .40 


O.ol 


A 

4 


1 

1 


n n^ i 

U.UO 


on n 
i yu.u 


i on n 
IZU.U 


n o 
u.z 


/I 7K 
4. / O 


q nQ 
o.Uo 


O OQ 

z.yo 





1 
1 


n c.n ~] 
U.OU J 


^n n 

L OU.U 


i on n 

IZU.U 


4.0 


n 

U.oO 


O 7Q 
Z. / O 


o 

Z.UO 


ft 


1 

1 


n i^n i 
U.OU J 


c;n n 
L OU.U 


i on n 
IzU.U 


n o 
U.z 


O QC 

z.oo 


n oq 
u.yo 


n O/i 
U.y4 


7 
I 


1 

1 


n c.n 
U.OU 


on n 
l yu.u 


i on n 
IZU.U 


A £ 
4.0 


1.00 


£2 OO 
o.Zz 


ft O/l 

o.y4 


o 
o 


1 

1 


U.OU 


on n 
yu.u 


i on n 

IZU.U 


n o 
u.z 


^ Q& 
O.oo 


o.oo 


o.oU 





o 
z 


n n^ i 

U.UO J 


c^n n 

L OU.U 


i on n 

IZU.U 


4.0 


n 

U.oo 


n c;7 
U.O / 


n Qft 
U.oO 


i n 

1U 




U.UO J 


c^n n 

L OU.U 


i on n 

IZU.U 


n o 
U.z 


n ft/i 

U.04 


n 07 
U.U / 


n 07 

U.U / 


ii 


2 


0.05 : 


L 90.0 


120.0 


4.5 


1.00 


3.04 


2.37 


12 


2 


0.05 ] 


I 90.0 


120.0 


0.2 


1.59 


0.34 


0.33 


13 


2 


0.50 : 


I 50.0 


120.0 


4.5 


0.49 


0.93 


0.69 


14 


2 


0.50 : 


I 50.0 


120.0 


0.2 


0.77 


0.10 


0.10 


15 


2 


0.50 : 


I 90.0 


120.0 


4.5 


1.21 


4.40 


3.71 


16 


2 


0.50 : 


I 90.0 


120.0 


0.2 


1.92 


0.49 


0.48 


17 


3 


0.05 ] 


L 50.0 


120.0 


4.5 


0.45 


0.78 


0.49 


18 


3 


0.05 ] 


L 50.0 


120.0 


0.2 


0.99 


0.17 


0.16 


19 


3 


0.05 ] 


L 90.0 


120.0 


4.5 


0.91 


2.50 


1.95 


20 


3 


0.05 ] 


I 90.0 


120.0 


0.2 


2.30 


0.71 


0.68 


21 


3 


0.50 : 


I 50.0 


120.0 


4.5 


0.47 


0.85 


0.63 


22 


3 


0.50 : 


I 50.0 


120.0 


0.2 


1.13 


0.22 


0.21 
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Table 4 — Continued 



Case 


Tvne a 




n b 




Ro 


t 


ly - L oelt 
(Ma,) 


Mn 
(10~ 5 M^ 

\ iU 1V± (BJ 


Mnn 

- Lv - L Uli. 


23 


3 


0.50 


1 


90.0 


120.0 


4.5 


0.99 


2.97 


2.51 


24 


3 


0.50 


! 


90.0 


120.0 


0.2 


2.68 


0.97 


0.94 


25 


4 


0.05 


! 


50.0 


120.0 


4.5 


0.23 


0.21 


0.13 


26 


4 


0.05 


! 


50.0 


120.0 


0.2 


0.52 


0.05 


0.04 


27 


4 


0.05 


! 


90.0 


120.0 


4.5 


0.44 


0.58 


0.45 


28 


4 


0.05 


! 


90.0 


120.0 


0.2 


1.20 


0.19 


0.19 


29 


4 


0.50 


! 


50.0 


120.0 


4.5 


0.24 


0.22 


0.16 


30 


4 


0.50 


! 


50.0 


120.0 


0.2 


0.59 


0.06 


0.06 


31 


4 


0.50 


! 


90.0 


120.0 


4.5 


0.46 


0.64 


0.54 


32 


4 


0.50 


! 


90.0 


120.0 


0.2 


1.40 


0.26 


0.26 


33 


1 


0.05 


2 


50.0 


120.0 


4.5 


0.89 


3.08 


1.95 


34 


1 


0.05 


2 


50.0 


120.0 


0.2 


2.34 


0.95 


0.88 




i 
i 


U.UO 





on n 


1ZU.U 






s no 


U.Ol 


36 


1 


0.05 


2 


90.0 


120.0 


0.2 


5.96 


4.76 


4.59 


37 


1 


0.50 


2 


50.0 


120.0 


4.5 


0.88 


2.98 


2.20 


38 


1 


0.50 


2 


50.0 


120.0 


0.2 


2.82 


1.38 


1.32 


39 


1 


0.50 


2 


90.0 


120.0 


4.5 


1.76 


9.32 


7.87 


40 


1 


0.50 


2 


90.0 


120.0 


0.2 


7.12 


6.80 


6.65 


41 


2 


0.05 


2 


50.0 


120.0 


4.5 


0.39 


0.58 


0.37 


42 


2 


0.05 


2 


50.0 


120.0 


0.2 


0.66 


0.07 


0.07 


43 


2 


0.05 


2 


90.0 


120.0 


4.5 


1.03 


3.21 


2.50 


44 


2 


0.05 


2 


90.0 


120.0 


0.2 


1.65 


0.36 


0.35 



-34- 



Table 4 — Continued 



Case 


Type a 




n b 






t 


M be lt 


M D 


M DR 














(Gyr) 


(M e ) 


(10- 5 M e ) 


(l(T 5 Me) 


45 


2 


0.50 


2 


50.0 


120.0 


4.5 


0.49 


0.95 


0.70 


46 


2 


0.50 


2 


50.0 


120.0 


0.2 


0.79 


0.11 


0.10 


47 


2 


0.50 


2 


90.0 


120.0 


4.5 


1.24 


4.68 


3.95 


48 


2 


0.50 


2 


90.0 


120.0 


0.2 


1.97 


0.52 


0.51 


49 


3 


0.05 


2 


50.0 


120.0 


4.5 


0.46 


0.82 


0.52 


50 


3 


0.05 


2 


50.0 


120.0 


0.2 


1.09 


0.21 


0.19 


51 


3 


0.05 


2 


90.0 


120.0 


4.5 


0.94 


2.68 


2.09 


52 


3 


0.05 


2 


90.0 


120.0 


0.2 


2.84 


1.08 


1.04 


53 


3 


0.50 


2 


50.0 


120.0 


4.5 


0.48 


0.90 


0.66 


54 


3 


0.50 


2 


50.0 


120.0 


0.2 


1.34 


0.31 


0.30 


55 


3 


0.50 


2 


90.0 


120.0 


4.5 


1.06 


3.37 


2.84 


56 


3 


0.50 


2 


90.0 


120.0 


0.2 


3.57 


1.71 


1.67 


57 


4 


0.05 


2 


50.0 


120.0 


4.5 


0.24 


0.22 


0.14 


58 


4 


0.05 


2 


50.0 


120.0 


0.2 


0.58 


0.06 


0.05 


59 


4 


0.05 


2 


90.0 


120.0 


4.5 


0.46 


0.63 


0.49 


60 


4 


0.05 


2 


90.0 


120.0 


0.2 


1.49 


0.30 


0.29 


61 


4 


0.50 


2 


50.0 


120.0 


4.5 


0.24 


0.23 


0.17 


62 


4 


0.50 


2 


50.0 


120.0 


0.2 


0.70 


0.08 


0.08 


63 


4 


0.50 


2 


90.0 


120.0 


4.5 


0.49 


0.72 


0.61 


64 


4 


0.50 


2 


90.0 


120.0 


0.2 


1.87 


0.47 


0.46 



a Table 1 constraint number. 
b See equation (3). 
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